Competition between superconductivity and nematic order in high-T c superconductor 
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We investigate the competition between superconductivity and nematic order in d-wave high-T c 
superconductor. Apart from the competitive interaction with superconducting order, the nematic 
order also couples strongly to gapless nodal quasiparticles. The interplay of these two interactions 
is analyzed by means of renormalization group method. An interesting consequence of ordering 
competition is the appearance of runaway behavior for some specific parameters, which implies an 
instability of first order transition. We show that the ratio between gap velocity and Fermi velocity, 
k = va/vf, of nodal quasiparticles plays a crucial role in determining the fixed point of the system. 
As k decreases, the possibility of first order transition is enhanced. At the nematic critical point 
where n — > driven by the critical nematic fluctuation, no stable fixed point exists and first order 
transition becomes inevitable. Our results indicate that gapless fermionic degrees of freedom should 
be taken into account in the theoretical description of competing orders. 

PACS numbers: 71.10.Hf, 73.43.Nq, 74.20. De 



I. INTRODUCTION 

Unconventional superconductors usually refer to the 
superconductors those can not be understood within the 
conventional Bardeen-Cooper-Schrieffer (BCS) theory. 
Notable examples of unconventional superconductors in- 
clude high-T c cuprate superconductor, heavy fermion su- 
perconductor, and iron-based superconductor. Unlike 
BCS superconductors, unconventional superconductivity 
is generally driven by electron-electron interactions, and 
often has a magnetic origin. Another interesting property 
of unconventional superconductor is that its ground state 
is not unique. In addition to the defining superconduct- 
ing state, unconventional superconductors also exhibit a 
variety of other symmetry-broken ground states, includ- 
ing antiferromagnetc, nematic, and stripe states, upon 
tuning such parameters as doping and pressure [l|45|. A 
widely recognized notion is that the long-range supercon- 
ducting order competes, and under certain circumstances 
coexists, with other long-range orders. The competition 
and possible coexistence between different orders can give 
rise to rich properties, and hence have attracted intense 
theoretical and experimental interest in the past years. 

The successful microscopic theory of competing orders 
has not yet been established to date, primarily because 
the pairing mechanism in most unconventional supercon- 
ductors is still undetermined. A realistic and commonly 
used strategy is to build low-energy effective field theory 
on phcnomcnological grounds. One can first write down 
the Ginzburg-Landau (GL) actions for two bosonic order 
parameters and then introduce certain coupling terms 
between these two scalar fields. Such generalized GL 
model has recently been applied to describe competing 
orders in a number of unconventional superconductors 
615]. An early success of such theoretical investigation 
is the prediction of field-induced antiferromagnetic core 
in the supercondcuting vortices of high-T c superconduc- 
tors [||. This prediction was subsequently confirmed in 
experiments [la, U3] ■ Interestingly, experiments further 



found that the antiferromagnetic order not only exists in 
the vortex cores, but also extends into the superconduct- 
ing region [l6| and exhibits nontrivial spatial modulation 
171 [l8| . A phenomenological field theory that contains a 



simple quadratic-quadratic coupling term between super- 
conducting and antiferromagnetic order parameters was 
put forward to understand these new findings @, Q ■ 

Recently, the issue of competing orders has attracted 
revived interest. It is found that the competitive inter- 
action between distinct orders can drive an instability, 
which gives rise to a general tendency of first order transi- 
tion (ill [Tij . This phenomenon may account for the first 
order transition observed in some unconventional super- 
conductors [![• In addition, nonuniform glassy electronic 
phases and Brazovskii type transitions are predicted to 
emerge due to competition between two long-range orders 
[To| . Another interesting observation is that the com- 
petition between superconducting and antiferromagnetic 
orders can help to judge the gap symmetry of iron-based 
superconductors [9|, ll2[. Furthermore, the competition 
between superconducting and nematic orders might be 
responsible for 15] the electronic anisotr opy observed in 
the vortex state of FeSe superconductor [19| . 

The effective field theory adopted in previous analysis 
of competing orders normally contains only two bosonic 
order parameters. The fermionic degrees of freedom are 
usually completely integrated out in the spirit of Hertz- 
Millis-Moriya (HMM) theory [2fj|-|22j. This integration 
procedure is expected to be applicable in systems that do 
not contain gapless fermionic excitations. For instance, 
the iron-based superconductors seem to have a s-wave 
energy gap, so the electronic excitations are fully gapped 
and can be safely integrated out 0, [l2j]. However, such 
integration manipulation is not always valid. Indeed, its 
validity has recently been questioned in several itinerant 
electron systems [23- 26] . In the systems that exhibit gap- 
less fermionic excitations, integrating out fermions may 
lead to singularities, especially in the vicinity of quan- 
tum critical point (QCP). Actually, infrared singularities 
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FIG. 1: Schematic phase diagram on (x, T)-plane of high-T c 
superconductors, x represents doping concentration. x\ and 
xi are QCPs of superconducting and nematic phase transi- 
tions, repectively. 



have been found on the border of several quantum phase 
transitions [23l - |26| . In order to properly describe the 
quantum critical behavior in these systems, it is more 
appropriate to maintain both bosonic order parameter 
and gapless fermions in the effective theory. When a 
long-range order competing with superconductivity also 
couples to gapless fermions, it would be interesting to 
go beyond HMM theory and examine the role of gap less 
fermions. Recent analysis presented in Refs. [l3|, Hjj did 
suggest nontrivial roles played by gapless fermions. 

In various superconductors, superconductivity com- 
pete with several possible orders. To examine the role of 
fermions, we wish to study a prototypical model which 
describes competition between two distinct long-range 
orders, contains gapless fermions, and in the meantime 
is technically controllable. In this paper, we choose to 
consider competition between superconductivity and ne- 
matic order in the contexts of high-T c superconductors. 
In recent years, there has been increasing experimental 
evidence pointing towards the existence of an electronic 
nematic phase in some high-T c superconductors 0-0, HB~ 
Hl|, especially YBa2Cti306+5 and E^S^CaC^Og+a. 
According to these experiments, a nematic order is pre- 
dicted to compete and coexist with superconductivity, 
which is schematically plotted in Fig. 1. The nematic 
phase transition and the coupling of nematic fluctuation 
with fcrmionic degrees of freedom have stimulated intense 
research efforts 0-0, H3, l32rj4l| . From a field-theoretic 
viewpoint, the nematic order parameter is a simple real 
scalar field and does not carry finite wave vector, which 
substantially simplifies theoretical calculations. 

It is known that high-T c superconductor has a d x 2_ y 2 
energy gap, which vanishes at four nodes, (±-|,±^). 
Therefore, gapless nodal quasiparticles (qps) are present 
even at the lowest energy in the superconducting phase. 
These nodal qps are believed to be responsible for many 
anomalous low-temperature properties of the supercon- 
ducting dome. When a nematic QCP exists somewhere 



in the superconducting dome, as shown in Fig. 1, the fluc- 
tuation of nematic order parameter will couple to these 
nodal qps. Such coupling can generate non-Fermi liquid 
behaviors and other unusual phenomena in the vicinity 
of nematic QCP [13, 136l - l4l| . It may also have significant 
effects on the interplay between superconductivity and 
nematic order, which is the topic of this paper. 

In the following, we first write down an effective field 
theory that describes both the competitive interaction 
between superconducting and nematic order parameters 
and the coupling between nematic order and nodal qps. 
We then perform a detailed renormalization group (RG) 
analysis |42| within such effective theory. In particular, 
we drive and solve the RG flow equations of all the phys- 
ical parameters in order to determine the possible stable 
fixed points. We show that the competitive interaction 
between superconducting and nematic order parameters 
still lead to first order transition after including gapless 
nodal qps. However, the nodal qps do have important 
impacts on the fixed-point structure. In particular, the 
velocity ratio k between gap velocity i>a and Fermi veloc- 
ity of nodal qps turns out to be a crucial controlling 
parameter. The region that exhibits runaway behavior 
is greatly enlarged as k decreases, implying an enhanced 
possibility of first order transition. At the nematic QCP 
where k — >■ driven by the critical nematic fluctuation, 
the system does not have any stable fixed point, so first 
order transition occurs inevitably. These results clearly 
indicate the importance of including gapless fermions in 
the effective field theory of competing orders. 

In Sec. UH we write down the effective action which 
contains two bosonic order parameters and gapless nodal 
qps. In Sec. IIII1 we make RG calculations and derive 
the flow equations for all parameters in the effective ac- 
tion. In Sec. IIV1 we present numerical solutions of the 
flow equations and discuss the physical implications. The 
paper is ended in Sec. fVl with summary and conclusion. 



II. EFFECTIVE FIELD THEORY OF 
COMPETING ORDERS 

We first need to write down an effective field theory to 
describe the competition between superconducting and 
nematic orders. This will be done largely on phenomeno- 
logical grounds. In the phase diagram presented in Fig.[TJ 
the horizonal axe is doping concentration x. The QCP 
of superconducting transition is x±, which is roughly 
x\ rs 0.05 in many high-T c superconductors. The antic- 
ipated QCP for nematic transition is represented by xi- 
So far, the precise value, and even the very existence, of 
X2 have not yet been unambiguously determined. Here, 
we assume that X2 is larger than xi, which implies a bulk 
coexistence of superconducting and nematic orders. 

In the present system, there are three types of degrees 
of freedom: superconducting order parameter ip, nematic 
order parameter 0, and gapless nodal qps The compe- 
tition between superconducting and nematic orders can 
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FIG. 2: (a): Free propagator of superconducting field ip; (b): 
Free propagator of nematic field <j>; (c): Free propagator of 
nodal qps <&. 



be described by a repulsive quadratic-quadratic coupling 
term, oc i/j 2 (f) 2 , which is widely adopted in the description 
of competing orders. In addition to this competitive in- 
teraction, the nematic order parameter <f> also interacts 
with gapless nodal qps <F, which is usually described by a 
Yukawa-type coupling term. There is, however, no direct 
coupling between the superconducting order parameter 
and nodal qps. First, the nodal qps are excited from the 
d x 2_ y 2 gap nodes where superconducting order parame- 
ter vanishes. Moreover, these qps are known to have a 
sharp peak and a very long lifetime in the superconduct- 
ing dome in the absence of competing orders [43j , so their 
coupling to ip must be quite weak. 

On the basis of the above qualitative analysis, we can 
write down the following partition function 



Z 



where the effective action is 
5 
Sij. 



5^ + 5^ + S-Jr + 5.00 

1 f d 3 q 

1 
2 
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(2tt) 3 
d 3 q 



(1) 
(2) 

d 2 rdr4> 4 , (3) 



(2tt) 3 v ' * ' T 2 



d 3 k 
(2^ 



[^[ii-iuj + v F k x r z + v A k v T x )^u 



+Vl i (-W + VpkyT* + V A k x T X )y 2 i], (5) 

5^0 = 7 J d 2 xdTip 2 (f) 2 , (6) 
5*0 = J d 2 xdT[\ cj ) ^\ l T x ^ ll + ^l l T x ^ 2i )], (7) 

where t x ' v,z are Pauli matrices and the flavor index i 
sums up 1 to N. <f{ represents nodal QPs excited from 
(§, § ) and (-§, -|) points, and the other two. The 
physical flavor of nodal qps, N — 2. Here, r is tuning 
parameter for nematic transition with r = at x 2 . vf,a 
are the Fermi velocity and gap velocity of nodal qps, 
respectively. The competitive interaction term ip 2 (f> 2 has 
a positive coefficient, 7 > 0. The propagators are shown 
in Fig. ([2|). In order to simplify calculations, it proves 
convenient to make two transformations |37| : <fr —> </>/Ao, 
and r — > Agr. It is now easy to rewrite Eq. ([7]) as 



5*0 = J d 2 xdT4>{^\ l T X ^ ll + ^l l T X ^ 2l )- 



(8) 



The effective action represented by Eq. (2) was stud- 
ied recently in Ref. 27]. It was demonstrated that both 



superfluid density and critical temperature T c are sig- 
nificantly suppressed at nematic QCP x 2 . However, 
the superconducting order parameter if) was assumed in 
Ref. [27j to be classical, which is valid only when x 2 is 
not close to x\. In this paper, we go beyond such approx- 
imation and consider the quantum fluctuations of both 
superconducting and nematic order parameters. 

As emphasized in Ref. [27|, the gapless nodal qps can 
have important impacts on the competition between su- 
perconductivity and nematic order. The simplest way to 
include the fermionic degrees of freedom is to introduce 
the polarization function H(q) due to nodal qps into the 
effective action of nematic order, S$. To the leading or- 
der of 1/iV-expansion, the polarization function H(q) is 
represented by the one-loop Feynman diagram shown in 
Fig. [3]^a) and formally given by [37], EH 

n(e,q) = 7V J ^J'ft[r B Go(w,k)r B Go(w + e,k + q)] ) 



where 



Go(w,k) 



f 



-iui + vpk x T z + v A k y T x 



is the free propagator for nodal QPs "Fi (free propagator 
for nodal QPs ^2 can be similarly written down). The 
polarization function n(e, q) has already been calculated 
previously [271 . |37| , and is known to have the form 



n(e,q) 



N 


e 2 + v 2 F q 2 x 




(e 2 + v 2 F q 2 + u|^)i/2 



(9) 



After including this term, the quadratic part of S$ be- 
comes 



(-2r + q 2 ) 



-2r 



n(<z)] 



(10) 



From the expression of polarization II(q), it is easy to 
see that inclusion of 11(g) does not change the dynamical 
exponent z = 1 of <f>. However, H(q) oc q, so it domi- 
nates over the kinetic term q 2 in the low energy regime. 
More importantly, the polarization Tl(q) introduces two 
important quantities, nodal qps' Fermi velocity vf and 
gap velocity va, into the effective action of (f>. 

Although the polarization H(q) represents the influ- 
ence of nodal qps, we can not completely integrate the 
nodal qps out and drop them from the effective theory. 
These gapless nodal qps should be maintained for sev- 
eral reasons. First, according to the general spirit of RG, 
one can safely integrate out high-energy modes at low 
energies. However, the nodal qps are gapless and hence 
exist even at the lowest energy. Integrating out gapless 
fermions completely may lead to unphysical singularities. 
Second, the coupling of gapless fermions with order pa- 
rameter fluctuation can lead to non- Fermi liquid behavior 
in observable quantities. In the specific case of nematic 
transition, the critical nematic fluctuation causes unusual 
fermion velocity renormalization and extreme anisotropy, 
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FIG. 3: (a): Polarization function for nematic field (f>; (b): 
Fermion self-energy correction due to nematic fluctuation. 



which would be overlooked if gapless fermions are fully 
integrated out. As will be shown below, the velocity ratio 
k = va/vf play nontrivial roles. 

The effective field theory contains seven parameters: 
a, r, /3, it, 7, vp,v/±. They are all subjected to renormal- 
izations due to the mutual interactions among three field 
operators: ip, <f), and ty. We will study the flow of 
these seven parameters under scaling transformations 
and eventually obtain seven RG equations. Since we 
study the t/j-cp interaction and ^-<f> interaction on equal 
footings, these seven RG equations are self-consistently 
coupled to each other. The low-energy behaviors of these 
parameters and the possible fixed points of the effective 
field theory can be determined by solving these coupled 
RG equations. 



III. RENORMALIZATION GROUP 
CALCULATIONS 

In this section, we make a RG analysis and obtain 
the flow equations of all the aforementioned parameters. 
In order to examine the impacts of gapless nodal qps, 
we go beyond the HMM theory and maintain nodal qps 
throughout our calculations. We first analyze the cou- 
pling between nematic order and nodal qps, and derive 
the RG equations for fermion velocities, vf,a- We then 
consider the competitive interaction between supercon- 
ducting and nematic order parameters, and obtain the 
RG equations for the rest five parameters, (a, r, /?, u, 7), 
which depend on the fermion velocities, vf,a- These 
seven equations are self-consistently coupled to each 
other since vp,A appearing in equations of (a, r, 0, u, 7) 
flow according to their own equations. 



A. Flow equations of vf and va 

The Yukawa-type interaction between nematic order 
and gapless nodal qps has been recently investigated in 
several papers [27l I36l - l4l| . It is well-known that the 
Fermi velocity of nodal qps is not e qual to the gap veloc- 
ity, i.e., vf 7^ va- Experiments (43l l44| have determined 
that the velocity ratio k = va/vf ~ 0.1. This ratio is a 
very important parameter since it enters into a number of 
observable quantities of high-T ff superconductors, includ- 
ing electric conductivit y 1451 l4q . thermal conductivity 
[46j |. superfluid density (46l |47|. and T c [47j. An interest- 



ing property revealed and discussed in Refs. [27J, |36| - |41| 
is that the velocity anisotropy is significantly enhanced 
by the nematic fluctuation. 

The calculation of nodal qps self-energy function and 
the derivation of flow equations have already been pre- 
sented in previous publications [13, EH , and therefore are 
not shown here. It is only necessary to summarize the 
basic calculations as well as the relevant results. To the 
leading order, the fermion self-energy is represented by 
the diagram Fig.[3Jb), and has the form 



ded 2 a 1 
E(w,k) = / -— ^Go(^ + e,k + q)- 



(2tt) 3 

As shown in Ref. [37|, it can be written as 

dS(k,o;) 
dlnA 

where 



n(«r 



(11) 



Ci(-iw) + C 2 v F k x T z + C 3 v A k y T x , (12) 



2tt 

d.r I deg(x,( 



c 2 = 



2(v A /v F ) 
Nir 3 

x 2 - cos 2 8 - {va/vf} 2 sin 
(x 2 + cos 2 9 + (va/vf) 2 sin 2 0) 2 

0/ / \ poo p2tz 

2{v a /vf) 



-00 JO 
2 „„c,2 a f„, /„, \2 „;„2 , 



dx / d0g(x,{ 
Jo 



Nir 3 

cos 2 9 — x 2 — (va/vf) 2 sin 2 • 
(x 2 + cos 2 9 + (va/vf) 2 sin 



2/n2 • 



C 3 



Q- 1 



2(va/v f ) 



Nir 3 
X 2 + cos 



dx 



d6G{x,( 





2 _l _2 a _ ( VA / VF y sin 2 Q 



(x 2 + cos 2 9 + (va/vf) 2 sin 



x 2 + cos 2 9 + (va/vf) 2 sin 
x 2 + sin 2 9 



x 2 + sin 2 6 + (va/vf) 2 cos 2 



2 a\2 ' 



(13) 



(14) 



(15) 



(16) 



Using Dyson equation, it is easy to get a renormalized 
fermion propagator 

G^^k.w) = -iio + v F k x T z +v A k y T x - S(k,w), (17) 

which leads to the following RG equations, 



dvF in n \ 

— = (Ci - C 2 )v F , 

dv A (n n s 

— = (Ci - C 3 )va, 



d(vA/vF) 

dl 



(C 2 -C 3 )(va/v f ), 



(18) 
(19) 
(20) 



where / > is running scale. A straightforward analysis 
showed that the ratio va/vf flows to zero at the low- 
est energy, giving rise to a novel fixed point of extreme 
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FIG. 4: One-loop corrections to mass parameters a in (a) and 
r in (b), respectively. 

velocity anisotropy (37j . Such fixed point in turn leads 
to a number of nontrivial consequences, such as unusual 
broadening of spectral function [36| , non- Fermi liquid be- 
havior [38[ , enhancement of dc thermal conductivity [39| , 
and suppression of superconductivity [27| ■ To analyze 
the influence of velocity renormalization and especially 
the extreme anisotropy manifested at ncmatic QCP on 
the nature of superconducting transition, we require that 
the constant fermion velocities appearing in the polariza- 
tion Tl(q) to flow with running scale I according to Eqs. 
CHI) and ©. 

The extreme velocity anisotropy is a special feature of 
the nematic QCP, where r = and the nematic fluctua- 
tion is critical. Away from nematic QCP, r ^ 0, so the 
nematic fluctuation leads only to relatively unimportant 
renormalization of fermion velocities. For r ^ 0, df,a 
and therefore their ratio va/vf remain finite. 

B. Flow equations of a, r, ft, u, and 7 

We next consider the competitive interaction between 
nematic and superconducting order parameters. Our 
analysis follows closely the scheme presented in a recent 
work of She et al. [Oj . It was assumed in Ref . [T3| that 
all the fermionic degrees of freedom can be integrated out 
and their effects can be represented by the dynamical ex- 
ponent z. Compared with Ref. [l4j . the main difference 
here is the inclusion of polarization Tl(q) in the effective 
action of nematic order (f>, which is supposed to reflect the 
influence of nodal qps. The corresponding (sub)action 
that describes ordering competition is 

Scorn = S.,f, + Sif, + S^tj), (21) 

where 

J m=l y ' 

s * = /!^ h2r+q2+f2+n(< ' )1 * 2 

+5/n%5^G>)'(5:«-)* 

J m=l K ' 



= 7 J 11 V(k 4 ,w J )0(q 4 ,e i ) 

xS 2 (ki + k 2 + qi + q 2 ) S (u>i + uj 2 + £i + £2) • 

Before performing a standard RG analysis within this 
action, it is convenient to rescale momenta and energy 
by A, i.e, k -)• k/A, u> -)• lj/A. 

According to the sprit of RG theory [3 H3, S3, SI El , 
we employ the following scaling transformations: 

h = k[e~ l , (22) 
w = Je~\ (23) 

^(k, W )=^(k', W ')e 5i/2 ; (24) 
q { = q[e-\ (25) 
e = e'er 1 , (26) 

0(q,e) = 0'(q',e')e 5Z/2 ; (27) 

where i — x,y. Applying these RG transformations, we 
will be able to derive the RG flow equations for the pa- 
rameters a, ft, r, u, and 7. 

Each field operator can be separated into slow mode 
and fast mode, i.e., 

ijj = ^J s +^f, (28) 
4> = 4>s+4>f- (29) 

After introducing an UV cutoff A, we can define the slow 
mode of superconducting order parameter as ip s — ip(k) 
with < k < e~ l A and the fast mode as if>f = iMk) 
with e~'A < k < A, using the formalism of Ref. |42j . 
Based on such modes separation, the effective action (|2"Tj) 
is decomposed into three parts: S s that contains only 
slow modes, that contains only fast modes, and S s f 
that contains both slow and fast modes. More concretely, 
we have 

Scorn = S -\- -(- S ^ 

= + S S ^ + S^) + (S$ + Si) + S*f : (30) 

where 

s * = 5/w Ha + k2 + " 2> * 2 

J m=l v ' 

J m=l v ' 
= 7 / 11 ?/'(k l ,w 4 )0(q l ,e l ) 

J i=l,2 ^ ' 

xS 2 (ki + k 2 + qi + q 2 ) 6 (oji + u 2 + e x + e 2 ) ; 
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of 



1 f d 2 kdoj 

2 J 

f T-r cPkmduj, 
+ 2 



-2a + k 2 + cj 2 ) ^/ 



5 



/ _ 



n 

m— 1 

1 /• d 2 qde 



-2r + q 2 + e 2 



11(g)] </> 2 



n 



d 2 q TO <ie 



m c2 



(2tt) 3 



S" 5 / = 



rf 2 kdw 3 , , , , t i i i \ 

^3 2 ^f^f^sVs + U<f>f<j>f(f)s<Ps) 



n 



d 2 k i du> l d 2 q^ l de i 2 , , 
_ <5 (k 1+ k 2 + q 1+ q 2 ) 

i=l,2 V ' 

x<5 (cji + w 2 + ei + £2) (jipsipsififfif 
+-fi>fi> f 4> s <ps + 4ryipfip s <f>f<f» B ) . 

After this decomposition, now the partition function can 
be rearranged in the following way, 



x exp 



(s^ + s; + s'^) + (sf + sf) + s-i 

J 1ty a V<f> a V*V* exp (5^ + 5; + S; 



J Vtfi.Vfa'DVV* exp (S; s + S£ + 5. 



(31) 



The next step is to integrate over all the fast modes, 
and obtain an effective action of slow modes. The func- 
tional integration can be performed using the standard 
diagrammatic techniques. The propagators for the su- 
perconducting order tp and the nematic order </>, shown 
in Fig. @, are 

G^u) = kz , L , — > (32) 



G^(q, e) 



1 



-2r 



•n(«)' 



(33) 



The polarization appearing in (q, e) reflects the pres- 
ence of gapless nodal qps. As already pointed out, 11(g) 
dominates over the kinetic term q 2 in the low-energy 
regime. In order to further simplify the nematic prop- 
agator, we consider the close vicinity of nematic QCP 
where r is very small. In this case, we are allowed to 
approximate the nematic propagator by 



G^(q, e) 



1 



-2r + U(q) 
1 



1%) (l - 



2r 



(a) 



(b) 



O 




FIG. 5: One-loop corrections to quartic coefficients (3 in (a) 
and u in (b), respectively. 



2r 



U(q) ip( q y 



(34) 



This expansion is apparently inappropriate when r is not 
small. For large r, it is better to consider H(q)/r as small 
expanding parameter, which will be discussed later. 

Now let us work in the spherical coordinates. We first 
define e = vq = vqcos9, qi = qsinOcostp, and g 2 — 
q sin 9 sin p, with v being the velocity of nematic order 
parameter <p. For convenience, we could extract constant 
v from the energy u of superconducting order parameter 
ip- Accordingly, the velocities of nodal qps, vf and va, 
can be divided by v. After this manipulation, vp and 
va become dimensionless. Now the polarization function 
can be written in the form, 



U{q,e,ip) = qD{B,ip), 
where the function 

1 



(35) 



D{6,p) 



16vpVA 



cos 2 9 + v 2 F sin 2 9 cos 2 1 



cos 2 9 + v F sin 2 9 cos 2 p + v\ sin 2 9 sin 2 p 



cos 



Vp sin 2 9 sin 2 p 



cos 2 9 + v 2 F sin 2 9 sin 2 p + v\ sin 2 9 cos 2 p 



is dimensionless. Before proceeding with the next calcu- 
lations, it is helpful to define 



/>7T /*27r 

F x = Fi(v f ,va) = / d6 / dp- 
Jo Jo (2tt 

F 2 = F 2 (vp,v A ) = d9 dip 



F 3 = F 3 (v F ,v A ) = d$ d v 



7T p2-K 



suit/ 


(2tt) 3 L>(£ 


K<p)' 


sin 9 




(2tt) 3 D 2 


0,<P) 


sin 6 




(2tt) 3 L> 3 


0,<P) 



(36) 
,(37) 
.(38) 



With these arrangements, we can now turn to calculate 
the one-loop contribution to the RG equations of all pa- 
rameters. 
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The diagrams contributing to a to leading order are 
shown in Fig. 0Ja). We perform the following calcula- 
tions, 



S[a] 



b d 3 



+ 7 



3 ip s (q)i/j s (-q) 



(2tt) 

1 J3„l 



(2^ 



6 d 3 



/. (2tt) 
9 



(27T) 



3 

^,( g )V,(-g)(-2o) 



-4 (1 + 2^+7(^+2^2) 



1 A - 3 V,(9)^(-9)e 5 ' /2 e 5 ' /2 (-2a) 



(2tt) 3 



x cxp 



2a" 



3/3 
2^ 



2 (1 + 2a) + 7(f\ + 2rF 2 ) 



d 3 q 
(2tt) : 



r^.(«)V.(-?)(-2o , ) ) 



(39) 



where 



a' = a cxp 2 Z 

' 2a 



3/3 

-^{l + 2a)+j{F 1 +2rF 2 ) 



Its derivative with respect to running scale I is 
3£ „ , „ x , 7 



da 



-^(l + 2a) + ^(F 1+ 2rF 2 ) 



(40) 



By calculating the diagrams shown in Fig. EJb), the 
flow equation for parameter r can be obtained similarly, 



dr 
~dl 



= 2r 



^ { F 1 +2rF 2 ) + -^(l + 2a) 



(41) 



In the absence of fermionic degrees of freedom, F\= F 2 = 
o^t, then the flow equation of a is identical to that of r 
[14| . Such an "exchange symmetry" is certainly broken 
by gapless nodal qps via the polarization II(<7) appearing 
in the effective action of nematic order d>. 




FIG. 6: One-loop corrections to coupling constant 7. 



It leads to 

dp 



9/3 2 



+ 2 7 2 (F 2 + 4rF 3 ) 



(42) 



By calculating diagrams shown in Fig. [Sfb), the RG 
equation for u is found to be 



du 
~dl 



r 7 2 



18m 2 (F 2 + 4rF 3 ) 



(43) 



Once again, the influence of gapless nodal qps is reflected 
in the functions F 2 3. 



E. 7 

Now we consider the flow of 7, which characterizes the 
strength of competitive interaction. The one-loop correc- 
tions to the term ji/j 2 ^ 2 has three diagrams, presented 
in Fig. (|()]). Following similar procedure, we eventually 
obtain 

- / n ^^£^(£^)^ 

J m=l V ' 

where 

7' = 7 expjZ- ^ +6u{F 2 + 4rF 3 ). 
+87 (Fi +F 2 2r + Fi2a) z}. 
The flow equation of 7 is therefore given by 



d 1 r r 3/3 



6u{F 2 + 4rF 3 ) 



-87 (Fx +F 2 2r + Fl 2a) }. 



(44) 



D. fi,u 

The one- loop corrections to j3 are depicted in Fig.[5ja). 
By paralleling the steps performed in Eq. (|39|) , we can 
similarly obtain 



s I ) 



IV. NUMERICAL RESULTS AND PHYSICAL 
IMPLICATIONS 

In order to specify the possible fixed point of the sys- 
tem under consideration, we need to solve the RG equa- 
tions. For later reference, it is useful to list all the RG 
equations obtained in the previous section, 



where 



/3' = /3cxp^l-- 



Q/3 2 

J^+^(F 2 +4rF 3 ) 



da 
dr 

Ti = 2r - 



^L(l + 2a) + l( Fl +2rF 2 ) 
^(F 1 + 2rF 2 ) + ^L(l + 2a) 



, (45) 
,(46) 
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d0_ 
~dl 

da 
II 

d<~) 



P- 



7T 2 

„,2 



27 2 (F 2 + 4rF 3 ) 



18u 2 (F 2 + 4rF 3 ) 



3/3 

+ 6u(F 2 + 4rF 3 ) 



d,VF 

~~dT 

dvA 
~df 



+87(^1 +F 2 2r + Fi 2a) } 
(Ci - C7 2 )«f, 

(Ci - c s )« A . 



(47) 
(48) 

(49) 
(50) 
(51) 



Compared with the case in which the action contains 
only two bosonic order parameters, the gapless nodal qps 
show their existence through the three functions -Fi ;23 , 
which are determined by fermion velocities, i>f,a- We 
now address how the velocity ratio affects the fixed point 
structure of the interacting system. 

As already pointed out, extreme velocity anisotropy 
with k — y is a special property of nematic QCP with 
r = 0. For finite r, the velocities receive only unimpor- 
tant renormalizations, and the ratio k remains finite. To 
make our analysis simpler, we first consider the case of 
finite r and solve the above RG equations by assuming 
two representative constant values of n. From the so- 
lutions, we will know qualitatively how the fixed point 
evolves as n decreases. We then consider nematic QCP 
with r = 0, and solve the flow equations (g5j |4"8"1 14"71 14"51 
l49l) self-consistently after taking into account the extreme 
anisotropy limit, k — > 0. 



A. va/vf = 0.1 

First, we consider a constant v\/vp = 0.1, which is an 
appropriate value for the velocity ratio in some realistic 
high-T c superconductors [H, Hj, |44| , In the presence 
case, the dimensionless velocity is set to vp = \ Q, and 
therefore ua — 0.1. Substituting them into Eq. (|5rJ|) . we 
have 



D(0,<p) = Z 



cos" 



sin 2 9 cos 2 tp 



and 



Fx 
F 2 
F 3 



dO 



d6 



dO 



2ti 



dtp 



suit 



2tv 



d<p 



2 71 



dip 



(27T 


*D(( 






sin 6 




(27T 


3 D 2 






sin 6 





{2ir) 3 D 3 {9,tp) 



0.053, (52) 
a 0.057, (53) 
a 0.062. (54) 




u, 



FIG. 7: Schematic diagram of flow trajectories in 
(j3,u) plane for bare velocity ratio va/vf = 0.1. 
The stable fixed point (SFP) is (a*, r* , /?*, it*, 7*) ~ 
(0.053547,0.043491,1.013065,0.757126,0.754569). When /3 
and/or u fall in the region left to and below the solid crit- 
ical line, the system under consideration exhibits runaway 
behavior, indicating the occurrence of first order transition. 



Now the RG equations (g5l | 
proximated as 



I S3 BH IH can be ap- 



da 
~dl 
dr 
~dl 
df3 
~dl 
du 

~dl 

d") 

~dl 



2a- 
2r- 

P- 
u- 
7 {l 

+87 (0.053 + 0.114r + 0.106a) | 



T 3 ^ (1 + 2a) + ^(0.053 + 0.114r) 

47T 2 2 

y (0.053 + 0.114r) + (1 + 2a) 



0/32 

Ar + 2 7 2 (0.057 + 0.248r) 

-y 2 

4r + 18u 2 (0.057 + 0.248r) 



3/3 



6u(0.057 + 0.248r) 



(55) 
(56) 
(57) 
(58) 

(59) 



8 V \J cos 2 9 + sin 9 2 cos 2 tp + 0.01 sin 9 2 sin 2 tp 
cos 2 9 + sin 2 9 sin 2 tp \ 
\J cos 2 9 + sin 2 9 sin 2 tp + 0.01 sin 2 9 cos 2 tp J 



By setting all these derivatives to zero, the fixed points 
could be determined. After numerical calculations, we 
find the following solutions (unphysical ones have been 
discarded) : 

1) (a* , r* , /3*, u* , 7*) 

(0.034935, o'.041791, 0.377768, 0.540608, 1.355735), 

2) (a* , r* , /?*, it*, 7*) s» 

(0.040467, o'.035673, 0.487222, 0.374853, 1.433849), 

3) (a* , r* , /?*, it*, 7*) s» 

(0.050708, 0.030099, 0.757159, 0.266773, 1.348298), 

4) (a* , r* , /3*, u* , 7*) « 

(0.053547, o'.043491, 1.013065, 0.757126, 0.754569). 

We need to find out the stable fixed 
point from these four solutions. By fixing 

a = a*,r — r*,7 = 7*, and making numer- 
ical analysis, we find that (a*, r* , (3*, u* , 7*) s» 
(0.053547, 0.043491, 1.013065, 0.757126, 0.754569) corre- 
sponds to the stable fixed point. In fact, we only need 
to address the physical properties in the vicinity of the 
stable fixed point. At a = a* , r — r* . 7 = 7*, it appears 
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that the flow trajectories of (3 and u show runaway 
behavior when the initial values of (3 and u are below 
the critical values f3 c pa 0.0835 and/or u c pa 0.0625, 
which are shown in Fig. ([7]). As extensively discussed 
previously [ll|i EH l4§ - i50| , such runaway behavior implies 
a first-order phase transition. 



B. va/vf = 0.01 

We then assume that va/vf — 0.01, and examine how 
the fixed point structure is modified as the ratio va/vf 
decreases. In the present case, we have 




25 



D(e,<p) = - 



cos 2 6 + sin 2 6 cos 2 ip 



4 V ^cos 2 9 + sin 2 9 cos 2 ip + 0.0001 sin 2 9 sin 1 



cos 2 9 + sin 2 8 sin 2 ip 



\/cos 2 6> + sin 2 9 sin 2 ^ + 0.0001 sin 2 cos 2 </? 



and then 
*i = 

^2 = 
^3 = 



sin# 



7r />2ir 

de I dip 

Jo 

d8 / dv? 
Jo 

^(2,)3^(^ ) 



(2tt) 3 D(^ 




sin 6 




(2tt) 3 I> 2 




sin 6 





0.005274, (60) 
a 0.000561,(61) 
a 0.000062.(62) 



Hence, the set of RG equations gSJ Hi HH 01 Hi) 
approximately becomes 



da 

HI 

dr 
~dl 
d£ 
dl 
du 
~dl 
dr, 



2a- 
2r- 

P- 



■^4 (l + 2a) + -(0.00527+0.00112r) 

47T 2 2 

— (0.00527+0.00112r) + -^r(l + 2a) 

2 47T 2 



■p^ 2 

7T 2 
'„ 2 



+ 2 7 2 (0.00056 + 0.00025r) 
18u 2 (0.00056 + 0.00025r) 



,(63) 
,(64) 
(65) 
(66) 



u r 



FIG. 8: Schematic diagram of flow trajecto- 
ries in (/3, u) plane for the bare velocity ratio 
=3Ja/vf = 0.01. The SFP is (a* , r* , /3* , u* , 7*) w 
^(0.053548, 0.437666, 1.013124, 76.671917, 7.590670). Com- 

) pared to the case of va/vf ~ 0.1 (see Fig. (7)), the region 
exhibiting runaway behavior is greatly amplified. It seems 
that a smaller ratio va/vf leads to enhanced possibility of 
first order transition. 



Among these points, only the first is stable. The cor- 
responding critical values of j3, u are j3 c ps 0.0833 and 
u c pa 6.32, which are shown in Fig. ©. 

Compared with the results obtained at va/vf = 0.1, 
we find that the probability of first-order transition 
rapidly rises as velocity ratio n decreases, because the 
region of runaway behavior is significantly enhanced. It 
is nature to speculate that first-order transition may be 
inevitable when the velocity ratio k vanishes. In order 
to judge this speculation, we finally consider the ncmatic 
QCP where r = 0. 



C. r = 

At nematic QCP with r = 0, the velocity is no longer 
finite: it becomes scale-dependent, and vanishes at the 
lowest energy. In this case, we should study the fixed 
point by self-consistently solving the following equations, 



■di=n 1 ' 



35 

4r + 6u(0.00056 + 0.00025r) 
7H 



+87 (0.00527 + 0.00112r + 0.01055a) }. (67) 



By applying the same methods used in the case of 
Va/vf — 0.1, we also find four solutions: 

1) (a*,r*,/3*,u*,7*) pa 

(0.053548, 0.437666, 1.013124, 76.671917, 7.590670), 

2) (a* , r* , p* , u*, 7*) pa 

(0.034916, 0.420692, 0.377357, 54.780339, 13.639865), 

3) (a* , r* , /3*, u*, 7*) pa 

(0.040451, o'.358983, 0.486613, 37.941187, 14.432612), 

4) (a* , r* , (3* , u* , 7*) pa 

(0.050767, o'.302544, 0.758586, 26.941770, 13.562402). 



da 
~dT 
d/3 
H 
du 
~dJ 
d"f 
~dJ 
dvF 

~dT 

dvA 
~df 



2a- 



3/3 



4n2 (l + 2a) + lF x 



= 0, 



"9/3 2 

7T 2 


4 2 7 2 ,F 2 


=0, 


2 


-18u 2 F 2 


-0, 



7 1 



nr2 



-6uF 2 + 87(^1 +Fx2a) 

{Cx -C 2 )v F =0, 
(Ci -C 3 )v A =0, 
which are obtained at the nematic QCP r = 



(68) 
(69) 
(70) 
= 0,(71) 
(72) 
(73) 
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FIG. 9: At nematic QCP with r = 0, the velocities «f,a and 
their ratio va/vf all flow to zero due to critical nematic fluc- 
tuation at the lowest energy. The system under consideration 
has no SFP in this case and exhibits runaway behavior in 
the whole plane of (/3,u), so first order transition becomes 
unavoidable. 



Since the equations of vf,va are relatively indepen- 
dent of other parameters, they could be first solved. 
It is easy to get v* F — 0, v A = 0, which implies that 
F\ = i<2 = -Fa = 0. The rest equations now become 



1 + 2a) = 0, (74) 

0, (75) 

1, (76) 
= 0, (77) 

0, p* = 



da 


= 2a 


3/3 




4tt 2 1 


(10 


= (3- 


9/3 2 


~dl 


" TT 2 


du 


= u - 




~dl 


9 


<h 




-5 


~di 


= -<{ 



7 = 
_ j_ 

22 • 



which have two solutions: 1) r* = u* = 
a* = 0; 2) r* = u* = 7* = 0, /?* = a* 

One can check that both of these two solutions are 
unstable. Hence, there is no stale fixed point at all at 
nematic QCP where the velocity ratio va/vf vanishes 
due to critical nematic fluctuation. The corresponding 
schematic flow diagram is shown in Fig. ([9]). The run- 
away behavior is present in the whole region spanned by 
(/3,u), indicating that first-order phase transition always 
takes place. 

Now let us discuss a little more about the role of nodal 
qps. Within the HMM theory, the fermionic degrees of 
freedom can be fully integrated out. If we do so in our 
case, then we would obtain an effective action that con- 
sists of solely two bosonic order parameters, ip and <f>. 
The influence of gapless nodal qps is only reflected by the 
polarization H(q) which contributes a term H(q)(j> 2 to the 
effective Lagrangian. The fermion velocities, vf,a, have 
to take their bare values and can not be renormalized, 
because the coupling between nematic order and nodal 
qps is neglected once the nodal qps are completely in- 
tegrated out. It is therefore not possible to incorporate 
the extreme velocity anisotropy driven by the critical ne- 
matic fluctuation into the theoretical analysis. However, 



as shown in the above calculations, such extreme velocity 
anisotropy does have significant effects on the RG trajec- 
tories (please compare Fig. 9 with Fig. 7). Apparently, an 
important indication of results is that gapless nodal qps 
do play an essential role and therefore should be taken 
into account in the effective theory of competing orders. 



D. Relatively large r 

In the above calculations, we have simplified the ef- 
fective propagator of nematic order parameter, G<f,(q), in 
terms of small quantity r/H(q). This expansion is appro- 
priate in the close vicinity of nematic QCP with small r. 
When r is relatively large, we consider Tl(q)/r as small 
parameter and expand G^(q) as 



G^(q, e) 



1 



-2r + n(?) 
1 



-2r(l-2M 
2r 



4r 2 



(78) 



We have performed RG calculations using this propaga- 
tor for relatively large r, and found that our qualitative 
conclusion reached above does not change. 



V. SUMMARY AND DISCUSSION 

In summary, we have performed a RG analysis within 
an effective field theory of the competition between su- 
perconductivity and nematic order in the context of d- 
wave high-T c superconductors. Our results confirm the 
qualitative conclusion of Refs. [3,[2l| that ordering com- 
petition can generically give rise to first order transition. 
Different from previous treatments, we go beyond the 
HMM framework and incorporate gapless nodal qps ex- 
plicitly in our calculations. It is found that the ratio k 
between gap velocity va and Fermi velocity vf of nodal 
qps plays a crucial important role in determining the RG 
flow diagram of the system. As k decreases from its bare 
value, the possibility of first order transition is enhanced. 
When k approaches zero at the nematic QCP due to the 
critical nematic fluctuation, first order transition become 
inevitable because the system exhibits runaway behavior 
in the whole space spanned by parameters u and /?, the 
quartic coefficients of superconducting and nematic or- 
der parameters respectively. These results indicates that 
gapless fermionic degrees of freedom can neither be ig- 
nored nor fully integrate out. 

The competition between superconductivity and ne- 
matic order is merely a very simple example of ordering 
competition. It is more interesting to study the competi- 
tion between superconductivity and antiferromagnetism, 
which is a fundamental issue in high-T c superconductor, 
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heavy fermion superconductor, and iron-based supercon- 
ductor. Compared with the case of competing nematic 
order, the interplay between superconductivity and anti- 
ferromagnetism is more complicated and exhibits richer 
properties. For instance, the antiferromagnetic order 
parameter is a complex scalar field and carries a finite 
wave vector Q which is often incommensurate. Further- 
more, the antiferromagnetic order parameter may acquire 
a nontrivial dynamical exponent, z ^ 1, due to coupling 
with gapless fermions, which would make RG calculations 
more involved [14| . Nevertheless, despite these techni- 
cal difficulties, the general formalism presented in the 
present paper can be applied to analyze the influence of 



fermionic degrees of freedom on the interplay between 
superconductivity and antiferromagnetism. 
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